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ABSTRACT 



Context. The 12.81 ^m [Ne n] line has recently gained interest as a potential tracer of gas in the tenuous surface layers of circumstellar 
disks and in outflow-related shocks. Evidence has been found for a proportionality between [Ne ii] emission and X-ray luminosity, 
supporting the hypothesis that X-rays are responsible for the required ionization and heating of the gas. Alternatively, ionization and 
heating by EUV photons and in J-type (dissociative) shocks has been proposed. 

Aims. The T Tau multiple system harbors three stars with circumstellar disks, at least one strong X-ray source (T Tau N), and regions of 
shocked gas in the immediate vicinity. ISO and Spitzer spectra revealed remarkably strong [Ne ii] emission, but because of insufficient 
spatial and spectral resolution those observations could neither pinpoint where in the system the [Ne ii] emission arises, nor identify 
the emission mechanism. We aim to clarify this by observing the system with enough resolution to spatially separate the various 
components and spectrally resolve the line emission. 

Methods. We performed high-resolution (i?=30000) spectroscopy of the TTau triplet at ~0f.'4 spatial resolution with VISIR at the 
VLT early February 2008. We spatially separated T Tau N from the southern close binary T Tau S, as well as the structures of shocked 
gas surrounding the stars. The individual southern components Sa and Sb remained spatially unresolved in our observations. 
Results. The dominant component of [Ne n] emission is centered on T Tau S and has a spatial extent of FWHM~ 1" 1 in a Gaussian fit. 
We detect spatially extended red-shifted emission NW of the system and fainter blue-shifted emission to the SE, which we associate 
with the N-S outflow from T Tau S. Only a small fraction of the [Ne ii] emission appears directly related to the X-ray bright northern 
component. Shocks may account for a substantial and possibly dominant fraction of the observed [Nen] emission. We estimate the 
total [Ne ii] flux to be 23±6 x 10"" Wm"^, in good agreement with the values measured by ISO in late 1997 and Spitzer in early 2004. 
Conclusions. Our observations show that outflows rather than the disk surface may dominate the observed [Ne n] emission in stars 
with strong outflow activity. We propose [Neii] emission in jets as a major factor causing the observed large scatter in the Lx vs. 



^[NcII] 



relation. We argue that T Tau S is the driving source of the T Tau "NW-blob" 
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1. Introduction 

Circumstellar disks are indispensable in the formation process 
of stars and are the birthplaces of planetary systems. Infrared 
observations of the refractory material ("dust") show that the in- 
gredients needed to form terrestrial planets and cores of giant 
gas planets are present and in place (see Natta et al. 2007, for an 
overview of dust in proto-planetary disks). The bulk of the disk 
mass is present in gaseous form, and is very difficult to observe 
due to its very low average opacity. 

Infrared emission in the [Neii] fine structure transition at 
12.81 /im has recently gained interest as a potential tracer of gas 
in the tenuous upper layers of circumstellar disks, or of small 
amounts of gas in debris disks. The intrinsically much weaker 
[Nem ] transition at 15.5 5 yum is observed only in exceptional 
cases ( Lahuis et al. 2007 1. If Ne^ is present and the gas is heated 
to >10^ K, the fine structure transitions are excited and we may 



observe [Ne ii] emission ( Glassgold et al. 2007 Pascucci et al. 
[2007) [Lahuis et al.|[2007) [Ercolano et al.||2008l l. Several candi- 
date mechanisms for the ionization and heating of the gas have 



* Based on observations performed at ESO's La Silla-Paranal obser- 
vatory under DDT programme 280.C-5035 



been proposed. Neon can be ionized via K-shell absorption of 
stellar X-rays, X-ray irradiation also heats the gas to several 
thousand Kelvin C Glassgold eraL][2007} [Ercolano et"aLl[2008) . 
Alternatively, EUV photons may ionize neon and absorption of 
radiation by small grains or PAHs heats the gas via the photo- 
electric effect ( Gorti & Hollenbach[ [2008 ). Strong, dissociative 
(J-type) shocks constitute a third possible ionization and heating 
mechanism (e.g. van den Ancker et aL|l999| ). 

[Glassgold et al.[ ( [2007| l model the efiisct of X-ray iiTadia- 
tion on the surface of a circumstellar disk, demonstrating that 
an X-ray source with a luminosity typical of young stars pro- 
vides significant ionization of Neon. The X-ray irradiation also 
heats the gas in the disk atmosphere to several thousand K out 
to a radius of ~20AU, beyond which there is a rather abrupt 
drop in gas temperature. Including mechanical heating from e.g. 
wind-disk interaction only makes a minor quantitative differ- 
ence. Therefore, the [Nen] emission is restricted to the inner 
-20 AU of the disk. 

To clarify the nature of [Ne ii] emission in YSOs there are 
currently two viable approaches. The first is a "statistical" ap- 
proach, in which [Nen] luminosities of a sample of sources are 
compared to other observables, such as the X-ray luminosity. 
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to search for coiTelations. Observations with the Spitzer Space 
Telescope have provided robust measurements of [Ne ii] fluxes 
of tens of young stars (e.g. Lahuis et al. 2007 ). Current space 
based observations are very sensitive, allowing relatively weak 
[Neil] lines to be detected, but leave the line emission spec- 
trally and spatially unresolved, limiting their use as diagnostics 
for the emission mechanism. The second approach is to perform 
detailed observations of individual objects, in which the [Ne ii] 
emission is spectrally or spatially resolved. Such observations 
are currently only possible with ground based instrumentation, 
strongly reducing sensitivity compared to space-based measure- 
ments, and are only feasible for comparatively [Ne ii] bright ob- 
jects. In this paper, we follow the second approach. 

In a pioneering study following the statistical approach, 
Pascucci et al. ( |2007| l suggested a correlation between the 



strength of the [Nen] emission and the X-ray luminosity in 
young stars surrounded by circumstellar disks. However, their 
sample constituted only 4 detections covering 0.2 dex in L[Neii] 
and Lx, leaving the proposed correlation tentative. Based on a 
somewhat larger sample, Espaillat et al. ( 2007] l cast doubt on 
the proposed relation between L[Nen] and Lx- They suggest that 
i'LNeii] may instead be correlated to the accretion rate and pro- 
pose EUV radiation from accretion shocks, rather than stellar 
X-rays, to be the main agent for irradiation of the disk surface 
and formation of the [Nen] line. In a study following the sta- 
tistical approach that we have conducted in parallel to the work 
presented here, we have extended the diagram to 33 objects with 
both [Ne ii] and X-ray detections, spanning ~2 dex in both L[Neii] 
and Lx (Giidel et al. 2009j in prep.), and do find a general 
trend of increasing L[Neii] with increasing Lx- However, large 
scatter is evident, casting doubt on the proposed direct relation- 
ship between [Ne ii] emission and X-rays, and arguing that the 
[Neil] generation in young stars is more complex. In particu- 
lar, sources with strong outflow activity are typically found to be 
over-luminous in the [Ne ii] line. 

Here we investigate the nature of the [Ne ii] in YSOs by per- 
forming a detailed study of the prototype of young, low mass 
stars: TTau ( iJoy||1945| [Ambartsumian||I947| [T949] l. The TTau 
system suits our objectives particularly well for several reasons: 
1) it is a multiple system with 3 known disk-bearing young stars 
within the central arcsecond, of which at least one is a strong X- 
ray source (TTauN); 2) it contains diffuse regions of shocked 
gas related to outflows on scales of several arcseconds; 3) it 
shows very strong [Ne ii] emission, allowing ground based stud- 
ies at high spatial and spectral resolution. 

Located in the Taurus molecular cloud at a distance of 148 pc 
( [Loinard et al.|200"7 l), the T Tau system consists of the optically 
visible northern component T Tau N and the "infrared compan- 
ion" TTauS approximately 0'.'7 to die sou th (Dyck et al.|1982 i. 
The latter is itself a close binary ( [Koresk o 2000) with a cur- 
rent projected separation between the components Sa and Sb of 
-a.'O or -19 AU (e.g. |Kohler et al.|2008| l. Mass estimates of the 
system members are ~2Mo, ~2.1 Mq and ~0.8Mo for TTauN, 
Sa, and Sb, respectively (e.g. Kohler et al.|2008 i. TTau suffers 
modest and time-variable line of sight extinction of Ay ~ 1 mag 
([Kuhi 1974; Mel'Niko v & Grank in 2005; Loinard et al. 2007). 
Towards Sb the extinction is estimated to be ~15 mag (jDuchenej 
|et al.|[2005) , the extinction towards Sa is substantially higher 
and may be several tens of magnitudes ( Duchene et al.|[2005| l. 
All three components show strong infrared excess emission, 
indicative of circumstellar disks ( |Mendoza V.||1966{ |Skemer| 
|et al.|[2008[ ). The heliocentric radial velocity of TTauN is 



-i-19kms"' ( [Hartmann et al.| 1986b, those of Sa an d Sb ai-e -i-22 
and H-21 kms ', respectively ( Duchene et al.|200 5V 

One of the three stars, TTauN, is a very bright, luminous, 
and variable X-ray source, the X-rays being predominantly of 
coronal origin. ( [Neuhauser et al.|I995[|Stelzer et al.|2000[[GudeI| 
et al. 2007]). The other two stars are probably also prominent X- 
ray emitters but the strong absorption by the high intervening 
gas columns makes their detection close to the bright TTauN 
very difficult (see Sect. |4TT]below; [Giidel et"aLp 007 ). The TTau 
system was also shown to exhibit strong [Neii] emission by 'van] 
|den Ancker et al.[ ([T999), who tentatively associate this emission 
with shocked gas in the NW blob. 

In this study, we seek to clarify the origin of the [Ne ii] emis- 
sion in the TTau system. All three known stars, as well as the 
extended regions of shocked gas, fall within a single spatial res- 
olution element of existing spectra made with ISO and Spitzer, 
leaving the origin of the emission and the underlying mechanism 
undetermined. We tackle the problem using high spatial (~0'.'4) 
and spectral (7?~30000) resolution ground-based observations, 
resolving the various components of the system. 



2. Observations and data reduction 

TTau was observed with the mid infrared imager and spec- 
trograph VISIR (Lagage et al. 2004), mounted on Melipal, 
the third of VLTs four 8.2 m Unit Telescopes. Our primary 
goal was to obtain high resolution spectra aro und the [Neii] 
CPi/i-^Pi/i) fine-Structure line at 12.81355 //m (Yamada et al. 
1985"). Spectroscopic observations were performed with 3 slit 
orientations, each executed during a separate night in February 
2008 and calibrated independently (see Table [TJ. Additionally, 
imaging in an approximately 0.23 fim wide filter centered on 
12.81 /im was performed, aimed at detecting bright [Neii] fea- 
tures at positions not covered by our spectroscopic observations. 
At a 3cr detection threshold of ~5xlO"'^ Wm"^ arcsec"^, we did 
not detect any such features at radii beyond 1" from the system. 
Our spectra are much more sensitive to [Neii] emission since 
the high dispersion strongly dilutes the vastly dominant telluric 
background and continuum dust emission. Note that the [Neii] 
emission contributes only ~5% to the system flux integrated over 
the spectral passband of our imaging filter, the remaining ~95% 
being continuum dust emission. In the remainder of this paper, 
we will restrict the discussion to the spectroscopic observations. 

In an accompanying paper (van Boekel et al. 2009, in prep.) 
we report on the imaging observations, announcing unexpect- 
edly fast (~0.25 mag increase in 4 days) continuum brightness 
variations in TTauS at 12.8 pm. We attribute these to fast vari- 
ations of the irradiation of the disk surface of TTauSa, due 
to variations in the accretion luminosity of the central source. 
Observationally, there is no evidence for a direct relation the be- 
tween the IR continuum variability and the [Ne ii] emission (see 
also Sect. [4311. 



2.1. HR spectroscopy around 12.81 fim 

Long-slit spectroscopy of the T Tau system was performed with 
a slit width of 0'.'4 and a spectral resolution of /?~30000. In or- 
der to test various hypotheses to the origin of the [Ne ii] emis- 
sion we observed the system with 3 different slit orientations, 
spatially covering the main components of interest (see Fig.[T]i. 
Slit 1 covers the northern and southern component, and is ori- 
ented roughly perpendicular to the Sa-Sb separation. Slit 2 cov- 
ers the northern component and the "NW blob", and incidentally 
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Fig. 1. Continuum subtracted spectra of the TTau system showing the 12.81 jum [Neii] emission line at various positions. The 
apertures over which the spectra were extracted are shaded grey and solid lines connect the individual spectra and apertures. 
Approximate radial velocities ar e determined by fittin g Voigt profiles, and are indicated in kms"' relative to the assumed radial 
velocity of T Tau Sa, (H-22km s \ Duchene et al. 2005 1. Note that the line profiles are generally much broader than the instrumental 
FWHM of a; lOkms-i. 



catches T Tau S as well, albeit with significant slit losses. Slit 3 
covers the southern component, and is oriented along the Sa-Sb 
separation. All slits cover some of the diffuse gaseous emission 
seen in deep AO-assisted near-infrared images (e.g. Herbst et al. 
[2007)[BecFet al.|2008t|Gustafsson et al.|2008| ). 

All observations were performed by ESO staff during the 
nights starting 2, 3, and 6 February 2008, and are summarized 
in table[T] A relatively large chop throw of 18" was used to 
avoid possible faint extended emission to be lost in the sky sub- 
traction procedure. Any sufficiently bright emission within 9" 
of the central sources would be detectable using our setup. A 
spectroscopic calibrator was observed along with each science 
observation for telluric and flux calibration. The atmospheric 
transmission was calculated using ATRAN ( Lord|[T992] l, a wa- 
ter column of ^2 mm was found to give a good match to the 
calibration observations (note that there is only 1 strong water 
line in the spectral range covered, which is right at the red edge 
of our spectra and does not interfere with the [Neii] emission). 
Other telluric lines in our spectra are due to O3 and CO2, a strong 
CO2 line at 12. 8 1224 yum is the main spectral feature interfering 
with the [Neil] line. Our telluric correction using ATRAN re- 
moves all features to the noise level of our calibration observa- 
tions (SNR a: 30). Taking the standard deviation in the total sys- 
tem response (including atmospheric transparency) determined 



from 6 calibration measurements, we estimate the absolute flux 
calibration of our spectra to be accurate to 5% (Icr). Radial ve- 
locity corrections have been applied to account for the Earth's 
orbital motion, and the indicated velocities for all spectra shown 
in this paper are with resp ect to the systemic ve locity of T Tau Sa 
(-1-22 km s"' heliocentric, 



Duchene et al. 



2005 1. 



VISIR HR spectra suffer from "fringing": a modulation of 
the flux with wavelength, with a peak to peak amplitude of about 
15%. We found this effect to be stable over periods of at least 
several months by comparing our February data with observa- 
tions obtained in June 2008. We fitted a spline profile to the cali- 
bration observations and divided the science observations by this 
curve. This procedure yielded satisfactory results, we estimate 
the amplitude of any possible remaining artifacts to be <2%. 



3. Results 

In Fig.|2]we show the spectra of the northern and southern com- 
ponents, extracted from slit 1 using the following method. At 
each wavelength, the profile along the spatial direction is as- 
sumed to be the sum of 2 Gaussians (one for TTauN and one 
for TTau S), of which the central position, width, and amplitude 
are fitted to best match the observations. The volume of both 
Gaussians yields the flux estimates for both sources. In terms 
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Table 1. Log of the spectroscopic observations with 
VLTA'ISIR. 
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of received energy, both spectra are dominated by continuum 
emission. The [Ne ii] line is spectrally resolved, and is clearly 
concentrated on the southern component. The extraction of 1- 
dimensional spectra from our data as shown in Fig. |2] provides 
a useful first impression of the [Neii] emission in the system, 
but discards much of the spatial information in our 2D long-slit 
spectra. 

The continuum radiation we receive is thermal dust emis- 
sion from the disks around TTau North and South. Its width 
in the spatial direction (FWHM of Gaussian fits) ranges from 
0'.'43 to 0'.'48 between the different measurements, compared to 
0'.'40 to 0'.'49 for the associated calibrators. We thus conclude 
that the continuum emission at 12.8 fj.m is essentially spatially 
unresolved for both the northern and southern component (see 
also Fig.[3]and Sect. |3.1| l. At 12.8/vm, the dust does not show 
spectral structure over the small wavelength range covered by 
our spectra (A/l~0.035 fim). Likewise, the spatial profile is con- 
stant over our small spectral range. This allows the continuum 
emission to be subtracted by fitting the profile in the spectral 
range where no [Neii] emission is detected, and thus to isolate 
the line emission as well as interpolate of the continuum spatial 
profile to the wavelengths where we see [Ne ii] emission. We use 
the continuum subtracted spectra in our analysis. 

3.1. Qualitative description of tiie observed [Neii] emission 

Contrary to the continuum dust emission, the observed [Nen] 
emission is clearly spatially extended. In Fig. [3] we show spatial 
profiles of the continuum dust emission and the [Ne ii] emission 
centered on T Tau S (after subtraction of the continuum|^as seen 
though slit 3, as well as the PSF obtained from the calibration 
measurement performed immediately after the science obser- 
vation. The continuum emission remains spatially unresolved, 
whereas the [Neii] emission is clearly spatially resolved. The 
spatial extent of the line emission is FWHM~ V.' 1 in a Gaussian 
fit in both slit 1 and slit 3 (N-S and NW-SE direction, respec- 
tively). 

While in Fig. [T] we show [Nen] spectra extracted in various 
apertures at key positions in the system. Fig. |4] shows position- 
velocity diagrams for each of the slits, displaying our data with 
continuous spatial sampling. The strongest [Neii] contribution 
arises in the I'.'l spatially extended component centered on 
TTauS. This emission is seen in apertures d, i, and h in Fig.[T] 
The E-W and NE-SW directions are not covered by our slits and 
thus the spatial extent in these directions cannot be determined. 
Note that the size of the emitting region, FWHM~l'.'l corre- 
sponding to FWHM~160 AU, is much larger than the disk size 
of either Sa or Sb: due to mutual tidal interactions, either disk 




-400 



-200 
V [km s"'] 



200 



Fig. 2. High resolution (Ra;30000) spectra of the north and 
south components of the T Tau system extracted from slit 1 (see 
Fig.[T]l. The spectra shown here were extracted using an approx- 
imate method (see Sect. [3]) and are useful to obtain a general 
picture of the continuum and [Ne n] line emission. For the more 
detailed analysis in this paper we use the continuum subtracted 
spectra shown in Fig. [T] The velocities are with respect to th e 



* In slit 3, we sampled the continuum blue-ward of the [Ne ii] line 
between -340 and -75kms"', and red-ward of the line between +100 
and -1-430 km s"'. 



radial velocity of TTauSa (h-22 kms , Duchene et al.||2005 i. 
The size of a spectral resolution element is plotted in the upper 
left corner. For reference, the atmospheric transmission is shown 
in the lower panel. Note that southern binary Sa-Sb remains spa- 
tially unresolved in these observations. 



cannot be larger than ~5 AU. The bright component centered on 
T Tau S has a velocity centroid that is blue-shifted compared to 
the stellar radial velocities. 

It is possible that the spatially unresolved disks of Sa and Sb 
contribute to the I'.'l extended component centered on TTauS. 
To estimate this contribution, we made a simple model of the 
line emission observed in slit 3, consisting of a point source 
and an extended, Gaussian component. Both components were 
convolved with the instrumental profile in the spatial direction, 
taken to be the profile of the unresolved continuum dust emis- 
sion, and compared to the observed profile. As expected, the 
extended component contributes dominantly to the total flux, 
and we estimate the possible contribution of the central point 
source(s) in TTau S to be <0.6 xlO"'^ Wm^^ 

Fainter extended emission is detected at larger distances, out 
to ~1'.'9 south of the TTauS and ~2'.'6 north of TTauN. The 
emission in the northern direction is systematically brighter than 
that in the southern direction. The extended emission seen in the 
S to SE direction (apertures e and j in Fig. [TJ is blue-shifted, in 
the N to NW direction we observe red-shifted emission (aper- 
tures b and g). We associate this emission with the N-S bipolar 
outflow discussed by Bohm & Solf ( I994| ), i.e. their "C" and 
"D" components. The red-shifted (*-i-40kms"') emission in the 
N direction comes to a halt (<-i-IOkms"') at the position of the 
"North- West blob" (apertures a and f). This is consistent with 
the existing notion, that the NW blob is a bowshock caused by 
an outflow impinging on ambient material. 

At the position of T Tau N (aperture c) we find a blue-shifted 
high-velocity component with a velocity of ^-125 km s"' . Both 
the location and the velocity of this emission match with that of 
the "B" component detected by Bohm & Solf ( 
associate this emission with a jet from TTauN. 



1994 1, and we 



Additionally, 
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a low-velocity blue-shifted component is seen at the position of 
T Tau N. Whether this emission can unambiguously be attributed 
to T Tau N is not clear, it may be part of the bright extended com- 
ponent centered on TTauS. Interestingly, the total [Neii] flux 
in this component (a; 10"'^ Wm"^) roughly equals the value pre- 
dicted by the tentative L[Neii] vs. Lx relation for TTau's very 
high X-ray luminosity of ^2x10"^' erg s"'. However, the central 
velocity of a!-22kms"'with respect to the stellar photosphere 
of T Tau N does not agree with a disk surface origin of this emis- 
sion. If directly related to TTauN, this emission may instead 
arise in a photo-evaporative flow from the T Tau N disk, of which 
we only see the approaching part since the receding, red-shifted 
part is obscured by the disk. 



Continuum and [Ne II] spatial profiles in slit 3 
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3.2. The origin of the N-S outflow and the "NW blob" 

Our [Ne ii] position-velocity diagrams (see Fig. [4| clearly show 
the bright emission centered on TTauS, as well as faint, ex- 
tended emission, particularly to the north. In the /9v-diagram of 
slit 1 we indicated the structures we associate with the north- 
ern part of the N-S outflow and the "NW blob". In both slit 1 
and slit 2 we can see that a connected structure is formed 
by the bright component centered on TTauS, material flow- 
ing northward with positive radial velocity, and the material in 
the NW blob that has approximately the stellar radial velocity. 
This strongly argues for TTauS being the source of the N-S 
outflow, as already proposed by Bohm & Solf] ( |I99 4) but later 



challenged by Herbst et al.| ( I996| l. It also shows that the out 



flow from T Tau S is indeed the driving force for the NW blob 
bowshock. 



3.3. Total [Neil] flux recovered 

The [Ne ii] line has previously been detected in the T Tau system 
using the ISO sateUite in late 1997 ( [van den Ancker et al.|1999| l, 
and with the Spitzer Space Telescope in early 2004. The mea- 
sured line fluxes were 28+7x10"'^ and 24±4xl0""' Wm"^, 
respectively, during these epochs. This suggests that the [Ne ii] 
emission is fairly stable on timescales of several years, though 
measurements at more epochs are required to confirm this. It is 
also conceivable that short term variations exist, induced by X- 
ray flares or shocks related to variable accretion in T Tau S. 

With our narrow slit, we do not cover the whole aperture 
through which the satellite spectra were taken, but we do sam- 
ple the key regions at least partially. Since our slit width of 0'.'4 
is much smaller than the measured width of FWHM=~l'.' I of 
the [Neil] emission centered on TTauS, this component suffers 
from large slit losses. Approximating its spatial intensity pro- 
file with a point symmetric Gaussian, we find that a total [Ne n] 
flux from this region of ~11 xlO"'^ Wm"^. If the actual spatial 
extent in the E-W direction is smaller than I'.'l, we will have 
over-estimated the contribution from this component. 

North- West of the triplet our slits cover an estimated one 
third of the region from which the [Ne ii] appears to arise, and we 
find a flux of ~8 xlQ-"' Wm^^ (this includes the "NW blob"). 
South-East of the triplet we estimate the extended emission to 
contribute ~2 xlO"'^ Wm"^, assuming our slits cover one third 
of the emitting region. About ~2.5xl0 Wm"^ of emission 
is seen at the position of T Tau N, of which about 40% stems 
from the high-velocity blue shifted component that is associated 
with a jet from TTauN, and the rest from the low-velocity blue 
shifted component. 



Fig. 3. Normalized spatial profiles of the continuum and [Ne ii] 
emission centered on TTauS, as seen through slit 3. The cuts 
were made at the approximate peak of the [Neii] emission at 
+ 10kms"' with respect to the Sa systemic velocity, and were 
integrated over 7 pixels along the dispersion direction (corre- 
sponding to a passband of 9.8xlO"'*yL(m or 23 km s"'). For refer- 
ence, the PSF profile extracted from a calibration observation is 
plotted with a dashed grey curve. 



In total, we find a flux of 23+6 xlO""" W nT^, where the rel- 
atively large error bar reflects the uncertainties induced by the 
correction for the incomplete spatial coverage of our data (for 
reference, the total line flux directly observed in the gray shaded 
apertures indicated in Fig.[l]is ~8 xlO"'* Wm"^). Thus, we ob- 
tain the same line flux as seen by ISO and Spitzer, and show that 
the [Ne ii] flux of the T Tau system was relatively stable over the 
1997-2008 period. 



4. Discussion 

Our main goal was to identify the spatial origin of the [Neii] 
emission in the T Tau system, with a possible relation between 
X-rays and [Ne ii] emission in the back of our minds (GlassgoU] 
|et al. 2007[|Pa"scucci et al.|2007| l. The first and most obvious con- 
clusion we can draw is that in the T Tau system there is no obvi- 
ous direct relation between [Ne ii[ emission and X-rays: whereas 
T Tau N is the dominant X-ray source, the [Ne n] emission arises 
partly in a bright concentration centered on T Tau S that may 
constitute the inner parts of the N-S outflow from this source, 
and partly in the more extended regions of this outflow as well 
as in its terminal shock, the NW-blob. 

Only a small fraction of the [Nen] emission, <10%, seems 
directly related to TTauN: ~1 xlO""' Wm"- is emitted in 
the spatially unresolved high velocity (-125 km s"') compo- 
nent which we associate with the jet of TTauN, another 
~1.5 xlO""' Wm"^ is seen at the location of TTauN but does 
not seem to be originating in the disk surface as judged by its 
radial velocity of -22 km s"' with respect to the photosphere of 
T Tau N. Instead, it may arise in a photo-evaporative wind from 
the disk of T Tau N, or even be an outskirt of the spatially ex- 
tended r.'l component centered on TTauS. 

Thus, while some of the neon emission seen in the TTau 
system may be X-ray induced, a direct relation between stel- 
lar X-rays and [Ne ii] emission from an irradiated disk surface 
( [Glassgold et al. 2007[ l clearly does not hold in the case of T Tau. 
The vast majority of [Nen] emission detected in the TTau sys- 
tem does not arise in the surface of any disk. 
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Fig. 4. Position velocity diagrams of the continuum subtracted spectra, showing the [Ne ii] emission in the T Tau system along our 
three slits. The positions of the continuum emission from TTauN and TTauS are indicated with dashed and dotted vertical lines, 
respectively. In the left-most diagram (slit 1) we indicated the features we attribute to the northern, red-shifted part of the north-south 
outflow, and the "NW blob". Contours are drawn at 1, 2, 3, 6, 9, 12, 15, 18, and 21 times the noise level. 



4.1. Mechanism responsible for [Nen] emission 

What mechanism may be responsible for the ionization and ex- 
citation of the neon atoms in the various parts of TTau sys- 
tem? Since the majority of the observed [Ne ii] emission arises 
in a spatially extended outflow, shocks constitute a primary can- 
didate. In typical jet densities of ~10^cm"^ the shock veloci- 
ties require d for produ cing substantial [Neii] emission are 70- 
100kms''(| Hollenbach & McKe e 1989 ). The line of sight veloc- 
ities we observe in the outflow from T Tau S are ~40 km s ' , but 
since the outflow is thought to be oriented close to the plane of 
the sky the actual velocities could plausibly be 100-300 km s"'. 
If the jet hits the inner "wall" of an outflow cavity under a graz- 
ing angle, only the velocity component perpendicular to the wall 
is relevant for the shock strength and the effective velocities are 
only a fraction of the jet velocity. All this considered, shocks pro- 
vide a plausible mechanism for ionizing the material we observe 
in the [Neii] line, but more detailed modeling of the emitting 
medium is required to be conclusive. The needed temperatures 
for excitation of the fine structure transition of several thousand 
Kelvin require shock speeds of 10-20 km s"' which are easily 
reached everywhere in the observed outflow. 

The second chief candidate mechanism for ionization and 
heating of the gaseous material is the absorption of stellar X- 
rays. The effectiveness of this mechanism is difficult to assess. 
Whereas the level and spectral shape of the X-ray emission from 
TTauN are relatively well known, virtually no X-rays are de- 
tected from TTauS. However, the high extinction toward both 
stars in the southern close binary allows for a substantial intrinsic 
X-ray luminosity of these objects. Sb suffers an estimated extinc- 
tion of Av~15 mag ( [Duchene et al.|2005] l; the extinction toward 
Sa is substantially higher, though its precise value is not well 
constrained. Adopting Ay = 15 mag and a conversion to the gas 
column density appropriate for standard int erstellar gas-to-dus t 
mass ratios, A^h = 2 x 10^' Ay cm^^mag"^ ( Vuonget al.|2003 i, 
we simulated the observed X-ray spectrum from T Tau S assum- 
ing that the intrinsic spectrum is identical to the one derived for 



T Tau N ( Giidel et al. 2007| l. We find that the bremsstrahlung por- 
tion of the spectrum peaks at 2-3 keV, while radiation at 1 keV 
is suppressed by 3 orders of magnitude. The total energy flux 
(in erg cm"^ s"') reaching Earth in the [0.3,10] keV interval 
is suppressed by a factor of 4.5. These suppression factors are 
applicable to possible X-ray emission from Sb, the attenuation 
of X-rays from Sa would be substantially stronger Given the 
approximate relation between stellar mass and X-ray luminos- 
ity for T Tauri stars (Telleschi et al. 2007| l, we expect Sb to 
be intrinsically less luminous in X-rays than T Tau N by a fac- 
tor of ~5. The observed X-ray flux from Sb would thus become 
~5% of that of TTauN, which agrees remarkably well with the 
6.5±2.7% that lGiidel eT ar(2007 ) derive for the faint diffuse ex- 
tention seen in the Chandra image at the approximate position 
of TTauS. Thus, it is quite possible that the southern stars have 
significant X-ray emission, and Sa could be as X-ray bright as 
T Tau N without being unequivocally detected in the combined 
X-ray spectrum of the N-S system (see e.g. [Giidel et al.||2007| 
who identify the main observed X-ray source with T Tau N). We 



note that the similar masses of T Tau N and Sa (e.g. Kohler et al. 
|2008[ ) indeed suggest similar X-ray luminosities. 

We have made an order of magnitude estimate of the level 
of [Nen] emission that would be induced by X-ray ionization 
in an extended volume of gas in the outflow from TTauS. We 
do not attempt to model the [Ne ii] emission in detail, but rather 
wish to establish whether or not stellar X-rays form a plausible 
ionization mechanism for the outflow material, and focus on the 
more extended emission at distances of >1". The source of X- 
rays could be TTauN or one of the southern stars. We assume 
the outflow gas to be hot enough for the 12.81 yum transition to 
be excited, and consider only X-ray ionization without addition- 
ally requiring heating by X-rays. We followed the methodology 
of Krolik & Kallman ( 1983[ l considering a stellar X-ray source 
irradiating a gaseous nebula, subject to continuum absorption 
and ionization. Ionization is to a large extent by photoelectrons 
produced by the primary photoionization. An average photoion- 
ization cross section was used for a gas with cosmic composi- 
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tion (2.7 X 10^22 cnr^ at 1 keV, |Igea & Glassgold||1999|l. For 



,-2 



the calculation of the emissivity of the [Ne ii] transition, we fol- 
lowed Glassgold et al. (2007 ) using recombination rates from 
|ShuU & van Steen berg (198^ and |Keady & Kilcrease| ( |2000| . 
We assuming a characteristic cylindrical outflow volume with a 
radial depth of 100 AU at a distance of 300-400 AU and a cross- 
section radius of 150 AU, and a gas density of 10^ cm""*, i.e. we 
mimic the volume of the T Tau NW blob. Furthermore, we take 
a bremsstrahlung X-ray spectrum approximating that of T Tau N 
(Lx = 1.5 X 10^' erg s ' in the 0.3-10 keV range, and a coro- 
nal electron temperature ofkT = l .5 keV, Giidel et al.|2007| l, an 
ambient outflow gas temperature of 5000 K, and an ambient ion- 
ization fraction of 0. 1 . We find luminosities of order 10^^ erg s"' 
in the [Nen] 12.8 fj.m line, corresponding to ~ 4 x 10"'^ Wm"^ 
of observed flux. In apertures a) and f) in Fig [T] we detect a 
total of 7.6 X 10"'^ Wm"^, which translates to approximately 
16 X 10"'^ Wm"^ after correction for the incomplete coverage 
of NW by our apertures. Thus, the line flux we predict is some- 
what less but within a factor of a few of what the observations 
show. Since we know neither the depth of the source nor the 
density, both of which are extremely important for accurate es- 
timates, our order of magnitude estimate is entirely satisfactory 
to show that X-ray absorption is a plausible ionization mecha- 
nism in the outflow. We note that this estimate is not based on 
self-consistent calculations; the production of ambient electrons 
(for the assumed ionization fraction of the gas) could be due to 
shock heating although the X-rays themselves help increase the 
ionization fraction as well. 

X-ray heating will be eff'ective relatively close to the stars; 
extrapolating the eflective outer radius of the [Nen] emis- 
sion in the |Glassgold et al.| ( |2007| l and [Ercolano et al.| ( |2008| l 
models of ~20AU to the X-ray luminosity of TTauN, which 
is approximately 10 times higher than that of the canonical 
models, we may expect X-ray heating to be effective within 
VTO X 20 AU < 70 AU. This simple extrapolation ignores opti- 
cal depth effects, due to the lower densities in the outflow com- 
pared to the disk atmosphere in the aforementioned models, X- 
ray heating may still be effective to somewhat larger distances. 
However, we detect [Neii] emission out to ~2'.'5 from TTauN 
in the north and south direction, i.e. ~370AU in projection. 
Heating by X-rays will not be effective at these distances, but 
as argued before shocks provide sufficient heating at every posi- 
tion in the outflow. 

In conclusion, both shocks and X-ray irradiation provide 
plausible means of ionizing the neon atoms in the outflow from 
TTauS. More detailed investigations are needed to establish 
which mechanism is dominant. Heating by shocks is probably 
very effective throughout the outflow, while X-rays are unlikely 
to contribute significantly to the heating at radii >100 AU from 
the X-ray source(s). A "hybrid" model, in which X-rays are re- 
sponsible for the ionization and shocks provide the heating re- 
quired to excite the fine-structure transition, may provide the 
most adequate description of the [Ne ii] emission in the outflow 
from T Tau S and those of young stars with strong outflow activ- 
ity in general. 



4.2. Variability of the [Nen] emission? 

The very fast continuum flux variations at 12.8/im detected in 
our VISIR data indicate a variable accretion rate in the central 
source on scales of a few stellar radii (van Boekel et al. 2009, 
in prep.), which is likely accompanied by variations in the out- 
flow rate. Such variations may plausibly lead to variations in the 



[Ne ii] emission of T Tau, which is outflow dominated, as we 
have shown here. However, as outlined in Sect. |3.3[ we find the 
[Neil] flux to be the same at epochs in 1997, 2004, and 2008, 
within the substantial uncertainties of the individual measure- 
ments. This suggests that the [Neii] flux is relatively constant 
over long periods. The continuum flux of TTauS, on the con- 
trary, is known to be strongly variable (e.g. |Ghez et al.||1991| l, 
and in fact it has varied by a factor of 2 between the three afore- 
mentioned epochs. Thus, there is no direct relationship between 
the observed [Ne ii] and IR continuum fluxes. Note that if vari- 
ations in accretion and outflow rate (as traced by the IR con- 
tinuum fluxes) would result in changes in the [Ne ii] flux, they 
would do so with a certain time lag, corresponding to the time 
needed for the out-flowing material to reach the positions where 
we observe the [Ne ii] emission. Given the relatively large spatial 
scales involved, this time lag would be ~1 year to tens of years, 
depending on the actual (de-projected) outflow velocity and the 
region under consideration. 

5. Summary and conclusions 

We have conducted high spatial (0'.'4) and spectral (/?=30000) 
resolution observations of the TTau multiple system, in order 
to reveal the origin of the strong [Ne ii] emission seen in space 
based observations in which the emission remained spatially and 
spectrally unresolved. The T Tau system contains 3 circumstel- 
lar disks, at least 1 strong X-ray source, and extended regions of 
shocked gas, thus providing several potential mechanisms pro- 
ducing [Ne ii] emission and an opportunity to test which is effec- 
tive. 

We find that the [Ne ii] emission is not concentrated on the X- 
ray bright northern component, but instead consists of a spatially 
extended (FWHM~l'.'l) concentration centered on TTauS that 
is likely the central part of the known N-S outflow from this 
source, and more diffuse emission associated with this same out- 
flow. The total [Neii] flux we derive from our February 2008 
observations is, within uncertainties, equal to the flux seen in 
ISO and Spitzer spectra taken in late 1997 and early 2004, re- 
spectively, suggesting that the [Ne ii] emission was not strongly 
variable over the last decade. 

Recently, Glassgold et al. ( 2007| l proposed that irradiation of 
the disk surface by stellar X-rays can lead to observable [Ne ii] 
emission (see also Ercolano et al. 2008 ). Indeed, this mechanism 
was shown to plausibly account for the [Ne ii] emission in sev- 
eral stars ( Herczeg et al.|2007 Pascucci et al.|20(T7 1. However, 
we here show that in the T Tau system, there is no direct rela- 
tion between the observed X-ray flux and [Neii] emission from 
a disk surface. Instead, the [Ne ii] emission arises in an outflow. 
While we could not establish the ionization and excitation mech- 
anism of the neon with certainty at every location in the system, 
we show that that both shocks and the absorption of stellar X- 
rays provide plausible ionization mechanisms for the outflow 
material. Shocks are the favored mechanism for providing the 
required heating of outflow material, though X-rays may con- 
tribute within ~I00 AU of the stars. 

Generalizing our result, we argue that young stars that ex- 
hibit outflows and shocks related to strong accretion activity 
can show strong [Ne ii] emission that is not directly related to 
the stellar X-ray emission in the fashion of the Glassgold et al. 
( 2QQ7\ model. We propose outflows to be an important [Neii] 
contributor, that may dominate over the disk surface emission 
at high outflow rates. This is indeed supported by our statisti- 
cal study of a large sample of T Tauri stars that show a correla- 
tion between the [Ne ii] luminosity and outflow parameters, and 
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also indicates that stars with jets reveal excessive levels of [Ne n] 
emission ( Giidel et al.|2009 in prep.). Lastly, we emphasize the 



importance of more spatially or spectrally resolved studies of the 
[Neil] line in young stars, such as that of Herczeg et al. (2007]l 
and this work, in order to search for signatures of outflows and 
Keplerian rotation in a disk surface. This is a challenging task 
since in the other stars in which the [Neii] line has been de- 
tected, it is typically fainter by at least an order of magnitude 
compared to the line flux observed in the T Tau system. 
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